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9 ABSTRACT 

10 Current variability of precipitation (P) and its response to surface temperature (T) are analysed using coupled 

11 (CMIP5) and atmosphere-only (AMIP5) climate model simulations and compared with observational estimates. 

12 There is striking agreement between Global Precipitation Climatology Project (GPCP) observed and AMIP5 

13 simulated P anomalies over land both globally and in the tropics suggesting that prescribed sea surface 

14 temperature and realistic radiative forcings are sufficient for simulating the interannual variability in continental 

15 P. Differences between the observed and simulated P variability over the ocean, originate primarily from the 

16 wet tropical regions, in particular the western Pacific, but are reduced slightly after 1995. All datasets show 

17 positive responses of P to T globally of around 2 %/K for simulations and 3-4 %/K in GPCP observations but 

18 model responses over the tropical oceans are around 3 times smaller than GPCP over the period 1988-2005. The 

19 observed anticorrelation between land and ocean P, linked with El Nino Southern Oscillation, is captured by the 

20 simulations. All data sets over the tropical ocean show a tendency for wet regions to become wetter and dry 

21 regions drier with warming. Over the wet region (>75% precipitation percentile), the precipitation response is 

22 ~13-15%/K for GPCP and ~5%/K for models while trends in P are 2.4%/decade for GPCP, 0.6% /decade for 

23 CMIP5 and 0.9%/decade for AMIP5 suggesting that models are underestimating the precipitation responses or a 

24 deficiency exists in the satellite datasets. 

25 

26 1. Introduction 
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27 The change in the global water cycle in a warming climate is a primary concern of society [ Meehl el al., 2007]. 

28 Model projections have indicated significant water cycle changes, with the intensification of extreme 

29 precipitation, the already wet areas getting wetter and the dry areas getting drier [Allan et al., 2010; Seager and 

30 Naik, 2011; Noake et al., 2012]. There is a robust physical basis for expecting precipitation (P) to increase in the 

31 global mean and in particular for regions of moisture convergence as surface temperature (T) rises, relating to 

32 energy and moisture balance constraints [Held and Soden, 2006; Mitchell et al., 1987; Midler and O ’Gorman, 

33 2011; Seager and Naik, 201 1]. Using multi-satellite observations, Liu and Allan [2012] assessed the consistency 

34 of the observed variability in P, and it was found that there is good agreement among data sets including GPCP 

35 (Global Precipitation Climatology Project) [Adler et al, 2008], SSM/I (Special Sensor Microwave Imager) 

36 [Wentz and Spencer, 1998; Vila et al, 2010], AMSRE (Advanced Microwave Scanning Radiometer - Earth 

37 Observing System) [Lobl, 2001], and TMI (Tropical Rainfall Measuring Mission (TRMM) Microwave Imager) 

38 over the tropical ocean and between GPCP and the TRMM 3B42 product [Huffman et al., 2007] over the 

39 tropical land (expected since both data sets use very similar gauge analyses and methodologies). Comparing 

40 climate model simulations with observations over the tropical oceans, Allan et al. [2010] found that the wet 

41 region (highest 30% of monthly precipitation values) is becoming wetter and the dry region (lowest 70% of 

42 monthly precipitation values) is becoming drier. However, results are sensitive to data sets and time period [Liu 

43 and Allan, 2012]. 

44 In the present study, we assess the current changes in global P simulated by historical scenarios from phase 5 

45 of the Coupled Model Intercomparison Project (CMIP5) and the atmosphere-only experiments (AMIP5) which 

46 are forced by realistic sea surface temperature (SST) and sea ice and radiative forcings. The aim of the present 

47 study is to evaluate how realistic and robust the models are in simulating the recent past, particularly over the 

48 satellite microwave measurement era. We assess the consistency and discrepancy between the simulations and 

49 the observations which has implications for the confidence in the projections of future climate change. 

50 

51 2. Data sets 
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52 We consider three observational data sets in the present study (GPCP, TMI and TRMM 3B42; Table 1). The 

53 GPCP is a global blended data set at 2.5° resolution containing land-based rain-gauges, sounding observations, 

54 microwave radiometers and infrared radiances [Adler et al., 2008]. The TMI data set only covers the tropical 

55 ocean from 40°N to 40°S at 0.25° resolution. The TRMM 3B42 covers the area from 50°N to 50°S at 0.25° 

56 resolution including both the land and ocean area but changes in ocean P are not considered realistic [Liu and 

57 Allan, 2012] because the existing AMSU-B algorithm failed to detect light rain over oceans, particularly in the 

58 subtropical highs [Huffman et al . 2007]; a corrected version is expected to be available soon. The data over the 

59 land region are consistent with GPCP observations. Observed T is the temperature at 2 m from the European 

60 Centre for Medium-range Weather Forecasts (ECMWF) INTERIM reanalysis [Dee et al., 2011] accumulated 

61 from six hourly 0.25° data interpolated from the original N128 reduced Gaussian grid (-0.7°). Blended T from 

62 the HadCRUT3 data set [Brohan et al., 2006] is also used for comparison purpose. Ocean (land) points are 

63 defined where all four neighbouring grid points are also ocean (land), aggregating from a high resolution 

64 (0.25x0.25 degree) land/sea mask; coastal grid points, which may be less reliable in the observational data (e.g. 

65 Huffman and Bolvin, 201 1), are excluded from the ocean-only and land-only comparisons in both models and 

66 observations. Details of the currently available CMIP5 historical experiments (12 models) and the AMIP5 

67 experiments (10 models) and their forcings are at http://cmip-pcmdi.llnl.gov/cmip5/. To ensure equal weighting 

68 from each model, we consider only one ensemble member from each CMIP5 and AMIP5 model to form 

69 composite CMIP5 and AMIP5 data sets (Table 1). 

70 3. Temperature and precipitation variations 

71 The deseasonalized T and P anomalies from ERA INTERIM, CMIP5, AMIP5 and satellite observations are 

72 plotted in Fig. 1. Mean P is also plotted in Fig. SI and listed in Table S2. The reference period is from 1988- 

73 2004 except for the TMI and TRMM data sets (1998-2004). Unlike the AMIP experiment which prescribes 

74 observed SST, the CMIP5 T simulations do not follow ERA INTERIM and have a large standard deviation 

75 since CMIP5 models generate their own ocean variability. The CMIP5 simulations contain realistic radiative 

76 forcings and can simulate cooling after the volcanic eruptions of El Chichon in 1982 and Mount Pinatubo in 
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77 1991 that are qualitatively consistent with AMIP5 simulations and observations (e.g. Fig. lc). The El Nino 

78 effect in 1988, 1998, 2005 and the La Nina effect in 1985, 1989, 2008 are clearly seen in the AMIP5 and ERA 

79 INTERIM T anomalies (Figs, lg-li). 

80 There is striking agreement between observed and AMIP5 simulated P anomalies over land both globally 

81 (Fig. le; r=0.6) and in the Tropics (30°N -30°S) (Fig. Ik; r=0.7). This suggests that prescribing the observed 

82 SST and realistic radiative forcings is sufficient for simulating interannual variability in land P. In general 

83 warmer years are associated with negative land P anomalies as noted previously [Adler el al., 2008; Gu el al., 

84 2007] and will be discussed in Section 4. 

85 GPCP displays greater P variation than both CMIP5 and AMIP5 globally (Fig. If) (the standard deviation of 

86 P from GPCP (-0.03 mm/day) is also higher than the individual models (-0.02 mm/day)) which is determined 

87 by the global and tropical oceans (Figs. Id and lj), though both AMIP5 and observations show positive phase 

88 correlations with T anomalies after 1995. To investigate the origin of these discrepancies, P anomaly 

89 differences between the AMIP5 ensemble mean and GPCP are calculated over the tropical ocean (Fig. 2a). The 

90 anomaly difference standard deviation (red line in Fig. 2a) is slightly reduced after 1995. 

91 Based on the periods of positive and negative area mean anomaly differences in Fig. 2a, the maps of mean 

92 anomaly differences are calculated for all positive (P + ) and negative (P ) AMIP5 minus GPCP anomaly 

93 composites over the period 1988-2008. The difference of P + -P is plotted in Fig. 2b. Regions of positive 

94 difference (the west and central south Pacific and western Indian Ocean) display a sign of variation that is 

95 consistent with the anomaly differences. This is further confirmed by plotting correlations between the local P 

96 anomaly difference time series and that of the tropical ocean mean (Fig. 2c). The regions that appear to 

97 contribute most strongly to the changes in AMIP5-GPCP anomaly differences are associated with the largest 

98 climatology difference between AMIP5 mean and GPCP P (Fig. 2d). 

99 There are a number of changes to the observed ocean data used in this study which may contribute to the 

100 discrepancy discussed above. For GPCP the switch from Outgoing Longwave Radiation (OLR) Precipitation 

101 Index (OPI) to Adjusted Geosynchronous Observational Environmental Satellite (GOES) Precipitation Index 
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102 (AGPI) in mid- 1987 is known to introduce an inhomogeneity in variance. The higher quality of the AGPI is the 

103 basis for examining changes starting in 1988 as well as 1979. Subsequent transitions between SSM/I sensors in 

104 1992 and 1995, and a change in aggregating the infrared data in 1996 are considered unlikely to provoke 

105 significant differences. As well, the GPCP shifts from low-orbit to geosynchronous-orbit IR data over the 

106 Indian Ocean in mid- 1998 (Huffman and Bolvin, 201 1). Removing the Indian Ocean (20°E-120°E) from the 

107 analysis improves the AMIP5-GPCP comparison much less than removing the West Pacific Ocean (Fig. S2b, 

108 S3b; Table S3), suggesting that the shift in Indian Ocean IR coverage does not introduce an inhomogeneity. 

109 Finally, the source of surface data used in the SST analysis shifts from Comprehensive Ocean- Atmosphere Data 
no Set (COADS) to Global Telecommunications System (GTS) in 1998 (Hurrell et al. 2008), reducing the surface 

111 data population available to provide calibration thereafter, but not obviously biasing the results. 

112 Natural changes may also influence the GPCP-AMIP time-series discrepancy. Both models and 

113 observational retrievals tend to exhibit different errors for different mean states of the atmosphere and therefore 

114 one might anticipate bias changes as the atmosphere changes. For example, the changing character of El Nino 

115 Southern Oscillation (ENSO) from an East Pacific (EP) to Central Pacific (CP)-dominated El Nino [Yeh et al., 

116 2009] may influence the statistical comparison of AMIP5 and GPCP since the climate simulation bias is 

117 strongest in the west Pacific. Indeed, the CP El Nino years (1990, 1994 and 2004) appear to correspond with 

118 negative AMIP5-GPCP in Fig. 2a. A related issue is the shift in the Pacific Decadal Oscillation in the mid- 

119 1990's. Changes in volcanic activity may also influence the GPCP-AMIP differences (large volcanic eruptions 

120 early in the record in 1982 and 1991) and this is another possibility to explore (e.g. Gu et al. 2007). Additional 

121 joint work by modelers and observationalists is needed to explicate the basis for the differences. 

122 Fig. 2e shows the scatter plot of the P anomalies between the AMIP5 mean and GPCP data sets over the 

123 tropical ocean, together with fitted lines (thick) over two periods (1988-1995 and 1996-2008). The correlation 

124 coefficient is -0.11 for 1988-1995 and is 0.23 over 1996-2008. The fitted lines between individual models and 

125 GPCP are also plotted in thin dashed line over these two periods: all models have positive and higher 

126 correlations over 1996-2008. The error source is quite complicated and merits further investigation but 
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127 nevertheless is suggestive of deficiency of the ocean observations prior to the introduction of the SSM/I F13 

128 data in 1995. It is expected that the comparison should be improved using the final version of GPCP 2.2 data 

129 [Huffman and Bolvin, 2011]. 

130 4. Precipitation response to surface temperature variation 

131 Precipitation response to the seasonal and interannual surface temperature variations are displayed in Figs. 

132 3a-3c and quantified in Tables 2 and SI. The relationships from CMIP5 and AMIP5 models are very close over 

133 the different regions analysed. For comparison purposes, unless stated otherwise, the data period used from now 

134 on is 1988-2005 for CMIP5, AMIP5 and GPCP data sets and from 1998-2008 for the TMI and TRMM 3B42 

135 data sets. 

136 The thick solid fitted lines denote statistically significant correlation (r) between P and T based on the two- 

137 tailed test using Pearson critical values at the level of 5% (dashed fitted lines denote correlations are not 

138 significant). The degree of freedom of the time series is calculated by first determining the time interval (t Q ) 

139 between effectively independent samples [Yang and Tung, 1998] but additionally assuming t G <12. (assuming 

140 that periods separated by 12 or more months are independent). 

141 Over the tropical ocean, the correlations between P and T are all positive. The precipitation change is ~3%/K 

142 for CMIP5 and AMIP5 simulations. It is 10%/K for GPCP P and ERA INTERIM T and 7.9 %/K if HadCRUT3 

143 T is used, close to 10.9%/K calculated by Adler et al. [2008] using an earlier version of GPCP. 

144 Negative correlations over the tropical land (-3.4 %/K for CMIP5 and -1.9%/K for AMIP5) are similar to 

145 GPCP (-3.1 %/K using ERA INTERIM T and -1.2%/K for HadCRUT3 T), but is smaller than TRMM 3B42 (- 

146 10 %/K for ERA INTERIM T and -1 1%/K for HadCRUT3 T) although this is for a short time period and most 

147 of the correlations are not statistically significant. Over the globe the GPCP dP/dT is positive and higher than 

148 the models (Table 2). 

149 The response over the tropical ocean and the tropical land is of opposite sign (Fig. 3d) for all datasets. The 

150 correlations are strong and significant (Table 2) and relate to ENSO [Gu et al., 2007], although monsoons must 
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151 also play a vital role [Hsu et al. 2010]. A similar relationship is also found between the global land and the 

152 global ocean (Fig. 3e). 

153 The strong relationship between GPCP and AMIP5 precipitation anomalies over the tropical land (Fig. 3f) is 

154 evident for the periods 1979-2008 (r=0.71), 1988-2008 (r=0.75) and 1998-2008 (r=0.74) but is weaker for 

155 AMIP5/TRMM 3B42 (r=0.35) over the 1998-2008 period. The agreement between the AMIP5 ensemble mean 

156 and GPCP data over tropical and global land is encouraging and suggests a strong control of ocean temperature 

157 on land precipitation as noted previously [ Gimeno et al., 2010]. 

158 5. Responses from wet and dry regions over the tropical ocean 

159 To further understand the source of discrepancy between tropical ocean P anomalies we now analyse the 

160 variability in terms of the monthly rainfall intensity distribution. Following Liu and Allan [2012], monthly 

161 precipitation is divided into percentile bins in ascending order of intensity and the anomaly time series of P 

162 averaged over the percentile bin is calculated. The anomaly time series of the area-weighted T over the tropical 

163 ocean is also calculated and the linear least square fit gradient, dP/dT, is computed. The percentage change 

164 (dP%/dT) is calculated by dividing dP/dT by the mean P for each bin over the reference period of 1988-2004. 

165 The dP%/dT and dP%/dt trend over the precipitation percentile bins are plotted in Figs 4a-4b and computed in 

166 Table 3. The non-linear scale of precipitation percentile is chosen since the higher percentiles contribute more 

167 to overall precipitation. The response is uncertain over the lower percentile bins, but is in general negative, 

168 consistent with Allan et al. [2010]. The wet region is characterized by positive dP%/dT in all data sets although 

169 the GPCP response is stronger. For dP%/dt, there is no physical reason to anticipate trends in tropical mean P 

170 unless there are associated trends in T or radiative forcings [Andrews et al., 2010]. The bin separating the 

171 positive and negative responses is around the 75% percentile for both calculations, consistent with previous 

172 analysis [Allan et al., 2010]. 

173 dP%/dT relationships over the wet (>75% precipitation percentile) region are positive and significant for all 

174 data sets. For GPCP data over the wet region, the change is 15%/K, around three times the model simulated 

175 responses and explains the discrepancy identified for the tropical ocean mean dP/dT discussed in Section 4. 
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176 Over the dry region the changes in P from models and GPCP data are quite consistent (~ -6%/K) when ERA 

177 INTERIM T is used (Fig. 4a). 

178 The precipitation anomaly time series over the wet and dry regions is plotted in Figs. 4c and 4d. The general 

179 trend is positive over the wet region and negative over the dry region despite the reduced trend in T since the 

180 1998 El Nino. The correlations between P over the wet and dry regions are -0.62 and -0.74 for CMIP5 and 

181 AMIP5 respectively and are significant. The GPCP variation in dry region P appears inconsistent with the 

182 AMIP5 ensemble after 1998 and is suggestive of a change in the sensitivity to light rainfall; the correlation 

183 between P over the wet and dry regions is insignificant (-0.12). For GPCP data, the precipitation trend over the 

184 wet region is 2.4 %/decade, close to previous estimates by Allan et al. [2010] but larger than CMIP5 and 

185 AMIP5 responses. Consistent with the tropical ocean mean comparison, correlation between GPCP and AMIP5 

186 P in the wet region is improved after 1995 (r=0.06 over 1988-1995; r=0.72 over 1996-2008). Conversely, over 

187 the dry regions of the tropical ocean, agreement between AMIP5 and GPCP data becomes poorer after 1995 

188 (r=0.38 over 1988-1995 and r=0.15 over 1996-2008). Over the dry region the CMIP5 and AMIP5 responses are 

189 substantially smaller in magnitude than GPCP but all data sets show a drying of the dry regions, though the 

190 correlations (r~0.3) are insignificant. 

191 6. Summary 

192 Current changes in precipitation over land and ocean are diagnosed from CMIP5 climate model simulations 

193 and compared with blended observations from GPCP and data from the TRMM satellite. Agreement between 

194 precipitation anomalies from GPCP and AMIP5 data set over the land (r~0.6) indicates that the atmosphere 

195 processes over the land are well represented by simulations including realistic SST and sea-ice changes and 

196 radiative forcings. Discrepancies between the observed and simulated tropical ocean P variability is found to 

197 originate primarily from the wet regions, in particular the west Pacific, but is reduced for the most recent period 

198 (1996-2008). However, differences over the dry regions of the tropical ocean are also evident and show poorer 

199 agreement between AMIP5 and GPCP data after 1995. This suggests that observed precipitation variability over 
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200 the ocean is sensitive to changes in the observing system; changes in ENSO character combined with model- 

201 satellite bias spatial signature may also influence the AMIP5-GPCP bias and trend differences. 

202 Despite the discrepancies, in all datasets considered, global and tropical ocean precipitation increases robustly 

203 with warming although observed responses appear stronger than those from models. Over the time period 1988- 

204 2005 the responses are 2.0%/K for CMIP5, 2.3 %/K for AMIP5 and 3-4 %/K for GPCP over the globe. Tropical 

205 ocean responses are larger but the responses over the tropical ocean and the tropical land are of opposite sign 

206 due to ENSO variability [ Gu et al., 2007]. There is a weak negative relationship between P and T over tropical 

207 land but the relationship between precipitation over the tropical land and the tropical ocean is strongly negative 

208 (r < -0.5). 

209 The analysis of precipitation change with temperature and with time show positive changes over the high 

210 precipitation percentile bins and negative change over the lower precipitation percentile bins, consistent with 

211 previous studies [Lau and Wu, 201 1]. Though the detailed precipitation changes still vary from model to 

212 observations and from model to model, the general characteristics of the precipitation variation and responses to 

213 the surface temperature variation are consistent. This supports the strong physical basis for expecting increased 

214 global precipitation with warmer surface temperatures due to energy constraints [Allen and Ingram , 2002] and 

215 for anticipating enhanced precipitation minus evaporation patterns due to moisture balance constraints [Held 

216 and Soden, 2006] and energy constraints [Muller and O ’ Gorman , 201 1]. However, further work is required to 

217 disentangle fast precipitation responses to radiative forcings from the thermodynamic responses [Andrews et al., 

218 2010; Ming et al., 2010; Wild et al., 2008] and to resolve the discrepancy between current interannual 

219 variability in observed and simulated tropical ocean precipitation. 
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Fig. 1. Temperature and precipitation anomaly time series relative to the reference period of 1988-2004 
over the global (a-f) and the tropical (30°S-30°N) (g-1) areas except for TMI and TRMM 3B42 from 
1998-2004. The black line is ERA INTERIM for temperature (a-c and g-i) and GPCP for precipitation 
(d-f and j-1). Shaded curves denote the CMIP5 and AMIP5 ensemble mean ± one standard deviation . 
Five month running means are applied. 

Fig. 2. (a) Time series of the area mean P anomaly difference (AMIP5 ensemble mean minus GPCP) 
over the tropical ocean, together with the five month running mean (thick black line) and the standard 
deviation over 1979-1995 and 1996-2008 periods (red), (b) the mean difference composite between 
positive anomaly months and negative anomaly months from 1988-2008 based on (a), (c) the correlation 
between the local anomaly difference time series and that from (a) over the period of 1988-2008, (d) the 
P climatology difference between AMIP5 ensemble mean and GPCP over 1988-2008 and (e) scatter plot 
of tropical ocean P anomalies between AMIP5 ensemble mean and GPCP over 1988-1995 and 1996- 
2008 periods , together with the fitted lines from AMIP5 ensemble mean and individual models. 

Fig. 3. Scatter plot of P and T anomalies (a-c) and P anomalies over the land and the ocean (d-e) from 
CMIP5/AMIP5 models and satellite -based observations and between AMIP5 and observed P anomalies 
over tropical land (f). Plotted linear fits are solid where significant at the 95% confidence level. 

Fig. 4. The change of tropical ocean precipitation with (a) tropical ocean mean temperature (dP%/dT) 
and (b) time (dP%/dt) over different precipitation percentile bins and precipitation time series over the 
wet (c) (>75% precipitation percentile) and dry (d) (<75% precipitation percentile) regions. Also 
displayed are CMIP5 and AMIP5 ensemble mean (solid line) ± one standard deviation (shaded area). 
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315 Solid symbols highlight significant correlations over the percentile bin and the time series is five month 

316 running mean. The seasonal cycle has been removed from all datasets. 

317 

318 
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Table 1. Data sets and their properties (rl is the first member of the model run). 


Data set 

Resolution 
Lat x Lon 

AMIP5 

1979-2008 

monthly 

CMIP5 

1979-2005 

Monthly 

BCC-CSM 

2.77° x 2.81° 


rl 

CanESM2 

2.77° x 2.81° 

rl 

rl 

CCSM4 

0.94° xl. 25° 


rl 

CNRM-CM5 

1.39° x 1.41° 

rl 

rl 

CSIRO-Mk3.6 

1.85° xl. 88° 


rl 

GISS-E2 

2.0° x2.5° 

rl 

r3 

HadGEM2 

1.25° xl.88° 

rl 

rl 

INMCM4 

1.5° x 2.0° 

rl 

rl 

IPSL-CM5A-LR 

1.89° x 3.75° 

rl 

rl 

MIROC5 

1.39° x 1.41° 

rl 

rl 

MPI-ESM-LR 

1.85° xl.88° 

rl 


MRI-CGCM3 

l.ll°x 1.13° 

rl 

rl 

NorESMl-M 

1.89° x 2.5° 

rl 

rl 

GPCP v2.2 

Combined observed precipitation from satellite and 

1979-2010 

rain gauges. Monthly data, global ocean and land, 
2.5° resolution. 

TMI v4 
1997 - present 

Monthly data, tropical ocean only (40°N -40°S), 
0.25° resolutions. 

TRMM 3B42 v6 

Tropical ocean and land (50°N -50°S), 0.25° 

1998 - present 

resolution. (TMI, SSM/I, AMSR), daily data. 

ERA INTERIM 
1979 - present 

6 hourly, global, 0.25 

0 resolution. 


HadCRUT3 
1979 -2011 

Monthly data, 5° resolution. 





325 

326 Table 2. Relationships of dP/dT and dPi and /dP ocea n over different region and time period. Significant correlation coefficient 

327 (r) at the 95% confidence level are marked in bold. Am is the error range of the gradient m. Values are in the round 

328 bracket when HadCRUT3 T is used and values in square bracket are the ranges of m from ensemble members. TMI ocean 

329 and TRMM 3B42 land datasets are combined for dPi an(] / dP ocean calculations. The values for each model runs are listed in 

330 Table SI. 

331 




dP/dT 

dPland /dPocean 

Data set 

Period 

Global 

Tropical ocean 

Tropical land 

Global 

Tropical 



m±Am 

(%/K) 

r 

m±Am 

(%/K) 

r 

m±Am 

(%/K) 

r 

m±Am 

r 

m±Am 

r 

GPCP v2.2 

1988-2005 

3. 8+0.5 


0.48 

10.3+1.0 

0.57 

-3. 1+0.9 

-0.22 

-0.36+0.08 

-0.30 

-0.81+0.09 

-0.52 



(3. 1+0.5) 


(0.43) 

(7.9+0.9) 

(0.52) 

(-1. 2+1.0) 

(-0.08) 





TMI ocean/ 

1998-2008 




15.5+1.5 

0.68 

-10.0+1.5 

-0.51 



-0.97+0.11 

-0.61 

TRMM 3B42 land 





(17.2+1.5) 

(0.71) 

(-11.1+1.9) 

(-0.46) 





CMIP5 

1988-2005 

2.0+0.04 


0.72 

3. 1+0.1 

0.51 

-3. 4+0.2 

-0.31 

-1.1+0.03 

-0.52 

-1.6+0.04 

-0.59 



[0.7 to 2.9] 



[1.4 to 4.4] 


[-13.4 to 0.6] 


[-2.2 to -0.23] 


[-3.2 to -0.4] 


AMIP5 

1988-2005 

2.3+0.06 


0.63 

3.0+0.17 

0.35 

-1.9+0.33 

-0.12 

-1.2+0.04 

-0.54 

-1.5+0.05 

-0.54 



[1.6 to 3.6] 



[-0.4 to 5.5] 


[-7.4 to 0.7] 


[-1.9 to -0.7] 


[-2.7 to -0.7] 



332 

333 

334 

335 

336 

337 Table 3. Tropical precipitation change with temperature and time. Correlation (r) is in bold when significant at the 95% 

338 confidence level. Values are in the round bracket when HadCRUT3 T is used and values in square bracket are the ranges 

339 of m from ensemble members. The T is the area mean over the tropical ocean (30°N -30°S ). 

340 


Data set 

Period 

dP„ et /dT 

dP dry /dT 

dP wet /dt 


dPjr/t 



m±Am 
(% /K) 

r 

m±Am 
(% /K) 

r 

m±Am 

(%/dec) 

r 

m±Am 

(%/dec) 

r 

GPCP v2.2 

1988-2005 

15+1.0 

(13±0.8) 

0.71 

(0.75) 

-5.9±3.1 
(-1 1.6±2.6) 

-0.13 

(-0.30) 

2.4+0.32 

0.44 

-3.5+0.76 

-0.30 

CMIP5 

1988-2005 

4.6±0.2 
[1.7 to 8.3] 

0.53 

-5.4+0. 3 
[-13 to 4] 

-0.31 

0.6+0.09 
[-0.1 to 1.3 ] 

0.42 

-0.5+0.12 
[-2.6 to 1.4] 

-0.29 

AMIP5 

1988-2005 

5.4±0.3 
[1.8 to 8.0] 

0.40 

-6. 0+0. 5 
[-15 to 1.4] 

-0.24 

0.9±0.16 
[0.3 to 1.6] 

0.36 

-1.5+0.28 
[-3.7 to 0.2] 

-0.35 


341 

342 

343 
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Month 


- INTERIM Ocean + Land 

INTERIM Ocean 


■ INTERIM Land 
Model Ocean + Land 
Model Ocean 
Model Land 


GPCP v2.2 

TMI 

TRMM 3B42 

CMIP5 BCC 

CMIP5 CanESM2 

CMIP5 CCSM4 

CMIP5 CNRM-CM5 
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AMIP5 MIROC5 
AMIP5 MPI-ESM-LR 
AMIP5 MRI-CGCM3-2H 
- AMIP5 NorESMI-M 
CMIP5 mean 
AMIP5 mean 


Fig. SI. Area weighted mean values of temperature (a-d) and precipitation (e-j) over the reference period of 1988- 
2004 except for TMI and TRMM 3B42 from 1998-2004. In the top two rows, the thick solid lines are from ERA 
INTERIM and the light colour patches are from CMIP5/AMIP5 all member mean ± one standard variation. For 
precipitation, the thick lines are from observations and ensemble means, the thin solid lines are from CMIP5 
datasets and the dashed lines are from AMIP5 datasets. 
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Fig. S2. Precipitation anomaly time series of AMIP5 and GPCP over the tropical ocean, (a) Re -plot of Fig. 
lj in the paper; (b) Without Indian Ocean (20°E-120°E removed); (c) Without Indian Ocean and Western 
Pacific (20°E-160°E removed). 
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Fig. S3. Time series of the area mean P anomaly difference (AMIP5 ensemble mean minus GPCP) over 
the tropical ocean, together with the five month running mean (thick black line) and the standard 
deviation of the anomaly difference over 1979-1995 and 1996-2008 periods (red), (a) Re -plot of Fig. 2a 
in the paper; (b) Without Indian Ocean (20°E-120°E removed); (c) Without Indian Ocean and Western Pacific 
(20°E-160°E removed). 
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Table SI. Relations of dP/dT and dP land /dP ocea n over different regions for each model run. Significant correlation 
coefficient (r) is in bold. Am is the error range of the gradient m. The data period is from 1988-2005. 


Dataset 

dP/dT 

dP land / dP ocean 

Global 

Tropical ocean 

Tropical land 

Global 

Tropical 

m±Am 
(% /K) 

r 

m±Am 
(% /K) 

r 

m±Am 
(% /K) 

r 

m±Am 

r 

m±Am 

r 


BCC-CSM 

1.8+0.13 

0.69 

1.7+0.45 

0.25 

-1.6+0.76 

-0.14 

-1.6+0.13 

- 0.65 

-1.8+0.16 

- 0.61 


CanESM2 

1.9+0.11 

0.76 

2.9+0.26 

0.60 

-6.3+0.67 

- 0.55 

-1.0+0.13 

- 0.48 

-1.8+0.20 

- 0.52 


CCSM4 

2.2+0.12 

0.79 

2.1+0.33 

0.40 

-1.6+0.69 

-0.16 

-1.3+0.12 

- 0.58 

-2.6+0.18 

- 0.71 


CNRM-CM5 

2.6+0.13 

0.81 

3.5+0.35 

0.58 

0.7+0.57 

0.08 

-0.18+0.09 

-0.13 

-0.78+0.11 

- 0.44 


CSIRO-Mk3.6 

0.7+0.13 

0.35 

4.2+0.27 

0.73 

-13.1+0.83 

- 0.73 

-2.2+0.09 

- 0.85 

-3.2+0.14 

- 0.85 

in 

Dh 

GISS-E2 

2.1+0.11 

0.80 

3.5+0.29 

0.65 

-2.4+0.43 

-0.36 

-0.62+0.10 

-0.39 

-0.9+0.10 

- 0.52 

§ 

HadGEM2 

2.3+0.11 

0.83 

1.3+0.36 

0.25 

-1.7+0.72 

-0.17 

-0.28+0.15 

-0.13 

-1.3+0.17 

-0.45 

u 

INMCM4 

2.2+0.25 

0.51 

3.9+0.38 

0.58 

-2.5+1.19 

-0.14 

-1.8+0.16 

- 0.62 

-1.7+0.19 

- 0.53 


IPSL-CM5A-LR 

2.2+0.11 

0.79 

3.0+0.44 

0.41 

-1.2+0.60 

-0.13 

-0.32+0.10 

-0.22 

-0.4+0.10 

-0.26 


MIROC5 

1.8+0.10 

0.77 

3.1+0.28 

0.61 

-3.5+0.57 

-0.39 

-1.1+0.09 

- 0.64 

-1.8+0.14 

- 0.67 


MRI-CGCM3 

1.8+0.18 

0.57 

2.0+0.73 

0.19 

-1.8+0.84 

-0.15 

-0.9+0.13 

- 0.45 

-1.4+0.13 

- 0.59 


NorESMl-M 

3.0+0.17 

0.77 

4.4+0.41 

0.60 

-2.5+0.52 

-0.31 

-0.8+0.10 

- 0.46 

-1. 2+0.1 1 

- 0.61 


CanAM4 

1.7+0.17 

0.57 

3.2+0.29 

0.61 

-7.5+0.99 

- 0.46 

-1.8+0.14 

- 0.65 

-2.8+0.30 

- 0.54 


CNRM-CM5 

2.5+0.19 

0.66 

3.2+0.50 

0.40 

-0.6+0.71 

-0.06 

-0.7+0.11 

-0.40 

-0.9+0.12 

- 0.46 


GISS-E2 

2.1+0.22 

0.54 

3.3+0.39 

0.50 

-3.2+0.80 

-0.26 

-0.6+0.13 

-0.32 

-1.0+0.16 

-0.41 


HadGEM2 

1.7+0.15 

0.61 

3.1+0.44 

0.43 

-4.2+0.94 

-0.29 

-1.5+0.11 

- 0.68 

-1.4+0.15 

- 0.55 

in 

Oh 

INMCM4 

2.4+0.16 

0.72 

2.8+0.54 

0.34 

0.9+1.00 

0.06 

-1.5+0.14 

- 0.60 

-1.8+0.17 

- 0.60 


1PSL-CM5A-LR 

2.8+0.24 

0.63 

4.5+0.61 

0.45 

-1.6+1.80 

-0.06 

-1.9+0.15 

- 0.66 

-2.5+0.18 

- 0.69 

< 

MIROC5 

2.1+0.18 

0.62 

2.5+0.40 

0.39 

0.5+0.80 

0.05 

-0.8+0.13 

-0.40 

-0.7+0.22 

-0.20 


MPI-ESM-LR 

1.9+0.18 

0.58 

2.0+0.70 

0.19 

-3.5+1.30 

-0.18 

-1.4+0.10 

- 0.70 

-1.7+0.12 

- 0.69 


MRI-CGCM3 

2.1+0.23 

0.54 

-0.3+0.54 

-0.04 

-1.5+0.85 

-0.12 

-0.8+0.17 

- 0.29 

-0.7+0.19 

- 0.25 


NorESMl-M 

3.7+0.18 

0.81 

5.5+0.75 

0.45 

-0.1+1.29 

-0.01 

-1.1+0.12 

- 0.53 

-1.6+0.13 

- 0.65 
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98 Table S2. Mean precipitation (mm/day) over different regions for observations and each model run. The mean period is 

99 from 1988-2004 except for TMI and TRMM 3B42 data sets from 1998-2004. 

100 



Dataset 

Global 

Tropical ocean 

Tropical land 

Tropical ocean 

Tropical ocean 






wet 

dry 

G 

o 

GPCP v2.2 

2.7 

3.0 

2.9 

6.8 

1.0 

c3 

> 

Sh 

D 

TMI ocean / 


2.7 

2.8 

8.4 

0.4 

C/2 

X> 

O 

TRMM 3B42 land 







BCC-CSM 

2.8 

3.6 

3.0 

8.3 

1.1 


CanESM2 

2.8 

3.7 

2.5 

8.6 

1.1 


CCSM4 

3.0 

3.7 

3.4 

8.3 

1.3 


CNRM-CM5 

3.1 

3.9 

2.9 

7.9 

1.8 


CSIRO-Mk3.6 

2.9 

3.7 

2.7 

9.2 

0.9 

in 

p_ 

GISS-E2 

3.2 

4.1 

3.5 

8.9 

1.5 

§ 

HadGEM2 

3.1 

3.9 

2.9 

9.5 

1.0 

U 

INMCM4 

3.2 

4.1 

3.2 

8.2 

1.9 


IPSL-CM5A-LR 

2.7 

3.5 

2.5 

7.9 

1.2 


MIROC5 

3.2 

3.8 

3.7 

8.2 

1.5 


MRI-CGCM3 

2.9 

3.6 

2.9 

8.7 

0.9 


NorESMl-M 

2.8 

3.5 

3.3 

8.1 

1.0 


CanAM4 

2.8 

3.7 

2.6 

8.5 

1.2 


CNRM-CM5 

3.0 

4.0 

2.8 

8.3 

1.7 


GISS-E2 

3.2 

4.2 

3.3 

9.2 

1.5 


HadGEM2 

3.1 

4.0 

2.8 

9.7 

1.1 

in 

Dh 

INMCM4 

3.2 

4.1 

3.6 

8.3 

1.9 

§ 

1PSL-CM5A-LR 

2.8 

3.4 

2.9 

7.4 

1.4 

< 

MIROC5 

3.2 

4.0 

3.7 

8.5 

1.7 


MPI-ESM-LR 

3.0 

3.8 

2.7 

9.1 

1.1 


MRI-CGCM3 

3.0 

3.7 

2.9 

9.3 

0.9 


NorESMl-M 

2.9 

3.6 

3.7 

8.0 

1.3 


106 Table S3. Correlation coefficients of the tropical ocean precipitation anomalies 

107 between GPCP and AMIP5 ensemble mean over different periods 

108 



1979-2008 

1979-1995 

1988-1995 

1996-2008 

All (Fig. lj) 

0.27 

0.03 

-0.20 

0.66 

No Indian Ocean 

0.42 

0.20 

0.17 

0.67 

No Indian Ocean 
and W Pacific 

0.80 

0.70 

0.70 

0.88 
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